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Abstract 

The Ferric uptake regulatory protein (Fur) is a transcriptional regulatory protein that functions to control gene transcription 
in response to iron in a number of pathogenic bacteria. In this study, we applied a label-free, quantitative and high- 
throughput analysis method, Interferometric Reflectance Imaging Sensor (IRIS), to rapidly characterize Fur-DNA interactions 
in vitro with predicted Fur binding sequences in the genome of Neisseria gonorrhoeae, the causative agent of the sexually 
transmitted disease gonorrhea. IRIS can easily be applied to examine multiple protein-protein, protein-nucleotide and 
nucleotide-nucleotide complexes simultaneously and demonstrated here that seventy percent of the predicted Fur boxes in 
promoter regions of iron-induced genes bound to Fur in vitro with a range of affinities as observed using this microarray 
screening technology. Combining binding data with mRNA expression levels in a gonococcal fur mutant strain allowed us to 
identify five new gonococcal genes under Fur-mediated direct regulation. 
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Introduction 

The Interferometric Reflectance Imaging Sensor (IRIS) is a 
photometric biosensiiig technology, which is designed on tlie 
principles of interferometry and has been applied as a microarray- 
based screening technique for studying the interactions of 
biomacromolecules [1,2,3,4]. IRIS can be broadly applied to 
characterize many types of biomolecular interactions, such as 
DNA-DNA hybridization, protein-protein binding, and protein- 
DNA interaction in a quantitative format [1,2,3,4]. Specifically, 
this technique has previously demonstrated utility for probing 
transcription factor interactions with arrayed oligonucleotides [5] . 
Bacterial transcription factors control tight regulation of gene 
expression in response to host specific environmental niches in a 
number of human pathogenic bacteria. One such pathogen is 
Neisseria gonorrhoeae, which is the causative agent of the sexually 
transmitted disease gonorrhea, one of the most common infectious 
diseases worldwide. 

Gonococcal colonization of mucosal surfaces requires tight 
regulation of gene expression in response to host specific 
environmental niches including iron-limited host environments. 
This adaptation is mainly achieved by the up-regulation of iron- 
acquisition components, which are repressed by a ubiquitous 



bacterial regulatory protein, the Ferric Uptake Regulator protein 
(Fur) [6,7,8]. Fur is a conserved 15-17 kDa protein comprised of 
an amino-termiiial DNA binding domain and a carboxyl-terminal 
dimerization domain and is present in a number of diverse 
bacterial pathogens [9,10,11,12]. In the presence of iron (II) or 
other divalent cations [9,10,1 1,13], Fur forms dimers and binds to 
dsDNA in a sequence specific manner to manipulate transcription 
initiation [9,10,11,13,14]. In its most basic form, the iron-bound 
Fur binds to the -10 and —35 motifs in the promoter region of a 
gene to prevent binding of RNA polymerase and leads to 
inhibition of transcription initiation [14]. For example. Fur 
represses transcription without an iron cofactor, a process termed 
fl/)o-Fur regulation that has been mainly demonstrated in 
Helicobacter pylori [15,16]. In its apo-iorm Fur has also been 
reported to function as an activator, auto-regulating the fur gene in 
Vibrio vulnificus [17]. The iron-bound Fur dimer has also been 
reported to function as a transcriptional activator via direct 
binding to defined promoter regions [8,18,19,20,21,22,23,24]. 
Moreover, Fur influences secondary regulatory components, such 
as small RNAs and histone-like nucleotide binding protein (H-NS), 
to regulate a subset of genes indirectly [25,26,27,28,29,30,31, 
32,33]. 
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The Fur regulon of jV. gonorrhoeae has not been well defined due 
to the lack of a mutant strain and an efiicient and rapid method 
to characterize the interactions between the gonococcal Fur 
protein and predicted Fur binding sequences. To define the 
gonococcal Fur regulon more efficiently, we developed a 
microarray-based, label-free quantitative method utilizing IRIS 
to screen sequence-specific interactions of Fur to dsDNA in vitro 
[5,34]. Combining this approach with bioinformatical analysis of 
the gonococcal genome and transcriptional analysis of JV. 
gonorrhoeae wild type, fur mutant and fur complemented strains, 
we identified new Fur regulated genes. \Vc also describe potential 
additional functional roles for Fur mediated regulation in jV. 
gonorrhoeae. 

Materials and Methods 

Prediction of Fur boxes 

The sequences upstream (—400 to +.50) of ATG, for each open 
reading frame (ORF) in the M. gonorrhoeae FA 1090 genome 
(Genebank No. AE004969.1) were retrieved using RSAT 
(http://rsat.ulb.ac.be/rsat/) [35] to represent the putative pro- 
moter regions of each gene. Seven experimentally determined Fur 
boxes in the two pathogenic Neisseria species, M. meningitidis and „¥. 
gonorrhoeae, were used as the patterns for Fur box predictions 
(Table SI). Matched sequences in the putative promoter regions 
were identified using fuzznuc module from EMBOSS by setting 
the mismatch rate as 0.4 [36]. Each Fur box in Table SI 
provided a list of matched sequences in the entire genome. 
Subsequently, overlapping sequences across the seven hsts were 
kept as the final output of predicted Fur boxes. 

Interferometric Reflectance Imaging Sensor (IRIS), 
substituting EMSA 

IRIS functions on the principles of interferometry [1,2,3,4,37]. 
Briefly, the optical path length of surface-reflected light is affected 
by changes in accumulated biomass on that surface. Diflerences in 
the optical path lengths (OPD) between the top layer and buried 
silica surface can be measured very precisely, which allows optical 
height information to be converted to accumulated mass (density) 
on tlu' surface". With this tc'chnifjuc, biomass accumulation can 
easily be quantified on a surface with a sensitivity of ~ 10 pg/ mm^ 
for a large field-of-view (FOV) to study multiple binding 
interactions simultaneously. 

Microarray preparation 

Amino-labeled forward DNA strands (50 |J,M) (Table S2) were 
mixed with excess complementary strands to create dsDNA 
through standard hybridization methods: mixed DNA solutions 
were heated to 80°C for 10 min and then allowed to cool to room 
temperature in a slow, controlled manner (<l°C/min). IRIS 
sensor chips were functionalized with a copolymeric surface 
chemistry to immobilize the amino-functionalized DNA strands in 
high density [38,39] as previously demonstrated. The dsDNA 
solutions were spotted onto the pre-functionalized sensor chips 
using a non-contact piezo-driven spotter, the sciFLEXARRAYER 
S3 (Scienion, Inc. Berlin, Germany). After a 12 hr immobilization 
period, spottc'd chips were washed with Saline Sodium Citrate 
(SSC) buffer in decreasing concentrations (2x, 0.2x, and 0. Ix) to 
produce a microarray of different covalently-attached dsDNA 
probes. IRIS images were taken for each chip to determine the 
initial dsDNA densities corresponding to each spot in the 
microarray. The chips were then incubated with 0, 200, 400, 
and 800 nM of purified gonococcal Fur protein [8,40,41,42] in 
2 mL of binding buffer (20 mM Tris-Cl pH 7.9, 5 mM MgClj, 



40 mM KCl, 0.125 mM MnClz, 2 mM DTT, 10% Glycerol, 
0. 19 |Xg/nL poly dIdC), respectively at RT for 2 hr. Subsequently, 
chips were washed three times (3 min for each wash) using 10 mL 
binding buffer followed by a 5 sec rinse in 0. 1 x binding buffer, 
and images were taken of the post-incubation array. 

Binding interaction analysis 

Increases in mass density due to Fur binding were determined 
by subtracting initial mass densities, determined for each spot from 
the pre-incubation images, from the final mass density measure- 
ments made after incubation [6] (Figure SI). On a spot by spot 
basis, the number of immobilized dsDNA molecules was 
determined from knowledge of the measured OPD, a known 
IRIS conversion factor for OPD (nm) to mass surface density (ng/ 
mm^) for dsDNA, and the molecular weight of each dsDNA 
oligomer. For the post-incubation measurements, an increase in 
the measured OPD was attributed solely to Fur binding and this 
difference was converted to the number of bound Fur dimers using 
simflar information as described above for calculating the number 
of immobilized dsDNA probes. In this way, the number of bound 
Fur dimers could be calculated for each spot within the microarray 
to effectively determine the number of protein dimers present per 
dsDNA molecule that was immobilized. For these measurements, 
between 6 and 10 spots were used to derive statistics on mass 
density calculations among the same condition. 

Purification of gonococcal Fur 

Purification of gonococcal Fur j)rot(;in was performed as 
described previously [8]. Briefly, the ORF of the fur gene was 
PCR amplified from genomic DNA of N. gonorrhoeae F62 and 
cloned into a PET 15b vector (Novagen, San Diego, CA) by the 
restriction sites Ndel and BamHI. The 6xHis-tagged Fur protein 
was over-expressed in a E. colifur mutant strain, HBMVl 19 [43] 
with 0.1 mM IPTG overnight at RT and purified using Ni- 
charged resin according to the manufacture's protocol (Novagen) 
[44,45] . Purified Fur protein was dialyzed against buffer [50 mM 
Tris-Cl, 500 mM NaCl, 100 \lM MnClg, 10% glycerol, pH 7.9]. 
The 6xHis tag was cleaved using a biotinylated thrombin 
(Novagen) and removed by flowing through another Ni-charged 
resin [40,46,47]. 

Electrophoretic mobility shift assay (EMSA) 

The probes used in EMSAs had the same sequences as those 
probes used in IRIS (Table S2). The dsDNA probes were 
obtained by mixing with its complementary strand at equal 
concentrations (5 |J,M), which were subsequentiy labeled with 
[Y-^2p]-ATP using T4 DNA kinase (Applied Biosystems, Ambion, 
Carlsbad, CA). The radio-labeled probes were purified using a 
G25 Sephadex QuickSpin column (GE healthcare, Pittsburgh, 
PA). EMSA was performed as described previously [8,41,42]. 
Labeled dsDNA Probes (12.5 nM) were incubated with purified 
gonococcal Fur in binding buffer (20 mM Tris-Cl pH 7.9, 5 mM 
MgCla, 40 mM KCl, 0.125 mM MnClj, 2 mM DTT, 10% 
Glycerol, 0.19 Hg/^L poly dIdC, 3.125 ng/jtL BSA) at RT for 
30 min. For cold competition assays, unlabeled competitor DNA, 
as indicated (50 to 1000 fold), was added to the reaction. Each 
reaction mixture was electrophoresed on a native 6% polyacryl- 
amide gel (acrylamide/bisacrylamide ratio, 37.5:1 [wt/wt]) (Bio- 
Rad, Hercules, CA), dried at 80°C on a filter paper for 1 hr and 
detected by autoradiography. The integrated density of bound and 
unbound DNA bands was quantified using ImageJ software and 
the percentage of bound DNA corresponding to a Fur concen- 
tration was calculated accordingly in each lane. The binding 
affinity (KD, equilibrium dissociation constant) of Fur to each 
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Table 1. Predicted Fur binding sequences within the promoter regions of W. gonorrhoeae iron-induced genes. ^ 



Predicted Fur box 


Gene 


Function 


Fur 

binding 


Fur-regulation^ 


TTA A A ATAT/"A A 1 1 1 A ATI" 
1 1 AAAA 1 A 1 uAA 1 1 1 AA 1 L 




Fe-S oxidoreductases family 1 


_4 


NE 


TATA I^Tf~(~(~ r A CCC A 1 1 1 1 
1 A 1 AA 1 LHjI^ALLIjA 1 1 1 1 




Fe-S oxidoreductases family 1 


_4 


NE 


TT A A A ATA/~A A f~ t" A "TT A Tf~ 

\ \ AAAA 1 AtjAALLA 1 1 A 1 C 


NGUOO/J 


phosphoglycolate phosphatase 


+ 


NR 


TATA A h K (~ K (~ t~ ATX/ — TT 

1 A 1 AAAAAA(jA(jLA I ICj I I 


r\(\ 1/^1 


cytochrome c4 


+ 


NR 


TATTATATA A A 1 1 1 1 1 A f~ 1" 

1 A 1 1 A 1 A 1 AAA 1 1 1 1 1 ACjL 


NGUOl DD 


hypothetical DNA binding protein 


+ 


Repressed 


(~ h fi A Af~ A Ai"/"ATA f — Mil/" 

(jAAAALAALuAI AL 1 MIL 


NGUUiUz 


hypothetical protein 




NE 


TT A A TT A A 1 — 1111/ — 1 1 1 1 A 

1 1 AA 1 1 AAC 1 1 1 1 u 1 1 1 1 A 


NGU0i04 


phenylalanyl-tRNA synthetase beta subunlt 


+ 


NR 


f~AAAAf~AAAA Af~f~f~ A 1 1 1 1 

LAAAALiAAAAACCCjA I I I I 


M^/^r\0 ~7"7 


NadC family protein, transporter 




NE 


TA/~A/ — TA A/ — TAT/ — TTATTT 
1 ALAL 1 AAu 1 A 1 L 1 1 A 1 1 1 


wi~ r\r\A "^c 


putative 3-methyl-2-oxobutanoate hydroxymethyltransferase 


+ 


Activated 


TATA ATA A A AT/" A ATT/ — TT 
1 A 1 AA 1 AAAA 1 LAA 1 1 L 1 1 


lNL]UUo4 1 


type III restriction/modification system modification methylase 


+ 


Repressed 


Ar"A ATA A Kt — III/ — 1 1 1 ATA 
ALAA 1 AAAu 1 1 1 L 1 1 1 A 1 A 


NLjUL)/ I 1 


alcohol dehydrogenase 


,4 
+ 


NE 


1 1 1 AAAAAAAA Te~ A A 1 1 1 1 

1 1 1 AAAAAAAA ILAAI 1 1 1 


NGDUoyy 


transcription elongation factor, GreA 


+ 


NR 


Tf~ AAAAAA/~AA Tf~f" AT A T(~ 

\ CdAAAAAAuAA I CCA 1 A 1 C 


NGUl loy 


chaperonin, HsIO 


+ 


NR 


T A A T A Te~ A A T A T A TT/~ A TT 

1 AA 1 A 1 CAA 1 A 1 A 1 1 CiA 1 1 


NGUl 2o4 


hypothetical protein 


+ 


NR 


ACAAAGAAGTATACTTCTT 


NG01419 


hypothetical protein 


+ 


NR 


TAATATAAGCGGCGGTATT 


NG01684 


hypothetical protein 




NE 


CCATACAACTATAI 1 1 1 1 1 


NG01738 


NADH dehydrogenase 1 subunlt M 


+ 


Activated 


TAACAGCACGCTCATTGTC 


NG01745 


NADH dehydrogenase 1 subunlt G 




NE 


TCAAATAAGAATCGTTATC 


NG01751^ 


NADH dehydrogenase 1 subunit A 




NE 


AAGAAAAAAGATGATTTTC 


NG01845 


305 ribosomal protein SI 2 


+ 


NR 


TATAAAAAGAAAGTATTTT 


NG01948 


hypothetical protein 


+ 


NR 


TAAAAGAAAACTCATTCTC 


NG01957 


putative export protein 


+ 


Repressed 



^Iron regulation as previously determined by microarray analyses of W. gonorrhoeae grown in defined medium CDM (-Fe) or CDM with 10 jiM ferric nitrate (+Fe) [50,51]. 
^Fur regulation was determined by comparing mRNA levels in N. gonorrhoeae WT, fur mutant and fur complemented strains grown under iron-replete and iron-deplete 
conditions at 1 h after addition of iron and desferal using qRT-PCR. NR; not regulated by Fur. Activated; mRNA level of the gene was decreased in the fur mutant strain 
under either iron-replete or iron-deplete conditions compared to that in the wild type strain. Repressed; mRNA level of the gene was increased in the fur mutant strain 
under either iron-replete or iron-deplete conditions compared to that in the wild type strain. NE, not examined; transcriptional regulation was not tested in this study. 
^Fur regulation of the genes as previously determined [8]. 

^In vitro Fur binding to the 500 bp upstream sequence of ATG of these genes was determined using EMSA and foot printing in previous studies [8]. 
doi:1 0.1 371 /journal.pone.0096832.t001 



DNA probe was calculated using GraphPad Prism Software by 
plotting the percentage of bound DNA (Y axis) with the Fur 
concentration pC axis) and using the equation Y= (Bmax • X)/ 
(KD -F X). In this equation, Bmax is the total number of receptors 
expressed in the same units as the 3^ values. 

RNA purification and quantitative RT-PCR 

jV. gonorrhoeae F62 wild type, fur mutant and complemented 
strains [48] were plated on GCB agar plates (Remel, Thermo 
Scientific) and grown overnight at 37°C with 5% COj and used to 
inoculate CDM medium [49] containing 0.042% NajCOa with 
shaking for 2 hr. The cultures were diluted into fresh CDM 
medium to an OD500 = 0. 1 and incubated for an additional 2 hr. 
Subsequently, 100 |J,M ferric nitrate or 150 |xM desferal (deferox- 
amine mesylate) was added to create iron-replete (-l-Fe) and iron- 
deplete (-Fe) conditions, respectively. Bacterial pellets were 
collected 1 hr after the addition of iron or desferal. RNA was 
purified using an RNeasy kit (Qiagen, Hilden, Germany) and 
treated with DNase I. Quantitative RT-PCR was carried out using 
a One-Step QuantiTect SYBR green RT-PCR kit (QIAGEN) on 
an ABI Prism 7700 sequence detection system (Applied Biosys- 
tems, Foster City, CA) and primers are listed in Table S4. A total 
of 12.5 ng RNA was used in each reaction. The relative mRNA 
levels of each gene were evaluated using the comparative cycle 



threshold (AACy) method (Applied Biosystems, User Bulletin Mo. 2: 
Relative Quantification of Gene Expression, 1997). The relative 
expression level of each gene was normalized to the endogenous 
16S rRNA gene and represented as the ratio to the expression 
level of the WT-I-Fe sample. The experiments were performed four 
times and results are presented as the mean ± the standard 
deviation. Statistics was performed using the Student's t test and a 
P<0.05 threshold was considered as significant. 

Results 

Prediction of Fur boxes in the genome of N. gonorrhoeae 
FA 1090 

Construction of the IRIS chips required a set of double strand 
DNA probes containing target sequences. To design the dsDNA 
probes, we first predicted Fur binding sequences (Fur boxes) 
throughout the genome using a bioinformatics approach. Neisseria 
Fur boxes previously defined by DNase I foot printing studies were 
utiUzed to avoid the mismatches caused by the species-specific 
nucleotides in Fur boxes of different bacteria (Table SI). A total 
of 198 Fur boxes in the intergenic regions throughout the 
gonococcal genome were predicted in this study (Table S5). Due 
to the fact that transcription of Fur regulated genes responds to 
iron availability, we focused on the iron-regulated genes identified 
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Figure 1. Identification of Fur binding to predicted Fur boxes by label-free IRIS screening. A set of —50 bp dsDNA probes containing the 
predicted Fur box in the middle of the probes were designed for IRIS screening (Table S4). Each dsDNA probe was immobilized on a chip to produce 
a spot with a diameter of approximately 100 |im. Three concentrations of Fur protein, 200 nM (red bars), 400 nM (blue bars) and 800 nM (orange 
bars), were incubated with the individually prepared arrays (in addition to a 0 nM control incubation (black bars) and the binding of Fur protein to 
dsDNA spots was measured. A mass increase for each spot was represented as a differential spot height (DSH). The known Fur boxes in fur, norB and 
nspA promoter regions were used as positive controls, and the aniA promoter region was used as negative control in this IRIS assay. (A) Differential 
spot height for each dsDNA probe after incubation with 0 nM, 200 nM, 400 nM and 800 nM Fur protein, respectively. (B) The number of Fur dimers 
bound per dsDNA molecule of each probe as calculated from initial and post-incubation mass density measurements. The gene designations of W. 
gonorrhoeae F62 were assigned according to their homologues in W. gonorrhoeae FA1090. 
doi:1 0.1 371 /journal.pone.0096832.g001 



via microarray analyses [50,5 1] of JV. gonorrhoeae cultured under 
iron-replete and iron-deplete conditions. Forty-five iron-repressed 
genes with predicted Fur boxes in their promoter regions (Table 
S3) were identified, which included 18 out of 21 of the previously 
identified Fur-repressed genes/operons in jV. gonorrhoeae [52]. This 
consistency validated our prediction of Fur binding sequences in 
the gonococcal genome. Twenty-two Fur binding sequences were 
predicted in the promoter regions of 2 1 iron-induced genes in this 
study (Table 1) including 2 previously identified Fur-activated 
genes, NGO07 1 1 and NGO 1751 [8]. 

Interactions of Fur with predicted consensus sequences 
in N. gonorrhoeae 

We next characterized in vitro binding of the Fur protein to 
predicted Fur boxes in the promoter regions of iron-induced genes 
using our high-throughput, label-free technology-IRIS. As shown 
in Figure lA, Fur binding to the probes containing a predicted 
Fur box demonstrated a concentration dependent increase in 
differential spot height (DSH). The Fur boxes in the promoter 
regions of fur, norB and nspA, which have different binding affmities 



to Fur (26.1 nM, 1.2 nM and 724.8 nM, respectively) were 
utilized as positive controls and the sequences in the aniA promoter 
as a negative control [8]. EMSA previously demonstrated that 
while the aniA gene contains a predicated Fur box, the Fur protein 
does not bind to this region [8]. Consistently, Fur binding to the 
known Fur boxes in fur and norB promoter regions reached 
4.8±0.5 and 5.4±0.9 nm increase in differential spot height, 
respectively, when 200 nM Fur protein was added (Figure lA). 
These mass density increases corresponded to 3.8 and 3.6 Fur 
dimers binding to one dsDNA molecule, respectively (Figure IB). 
The differential spot height of the nspA Fur box probe continued to 
increase with Fur protein concentrations from 2.8±0.4 nm, 
3.5±0.1 nm to 4.4±0.1 nm (Figure lA). Two Fur dimers, 2.6 
dimers and 3.0 dimers/ n^pA Fur box were detected when 200 nM, 
400 nM and 800 nM Fur protein were added to the chips, 
respectively (Figure IB). As the negative control, the aniA 
promoter sequences showed a small mass density increase for 
the three concentrations of Fur protein (DSH< 1 nm) (Figure lA). 
Therefore, we defined a threshold value for a DSH of & 1 nm as a 
cutoff for determining specific Fur binding to dsDNA probes. 
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B NGO0073 NGO0101 




Kd=39.8±6.3 nM Kd = 406.1 ±144.8 nM 



rmp-negative control fur-positive control 




Ko=31.8±7.9 nM 



Figure 2. Fur binding to the predicted Fur boxes. (A) Cold competition assay for the specificity of Fur binding to the predicted Fur 
boxes. ^^P-labeled —50 bp dsDNAs were analyzed after incubation with Fur and unlabeled dsDNA probes (cold probes). Two types of cold 
competitor probes were used: fur, which contains a Fur box, and rmp, which does not bind to Fur. When the unlabeled probes fur compete out 
binding of labeled probes, but the unlabeled rmp cannot, the binding of Fur to the predicted Fur box was considered specific. Lane 1, free ^'^P-labeled 
DNA; Lane 2 through Lane 10 contained gonococcal Fur. For the fur and NGO0073 Fur box, Fur protein was added at a concentration of 100 nM and 
for the NGO0101 Fur box, 400 nM of Fur protein was added. The fold excess of the cold probes was increased from Lane 3 to Lane 6 and Lane 7 to 
Lane 1 0 ranging from 50 fold, 500 fold to 1 000 fold (indicated by the triangles). (B) Fur binding affinities determined by EiVISA. The ^^P labeled 
dsDNA probes were incubated with a gradient of concentrations of purified gonococcal Fur protein (Lane 1 to Lane 10, 0 nM, 5 nM, 25 nM, 50 nM, 
100 nM, 200 nM, 400 nM, 600 nM, 800 nM and 1000 nM, respectively). Arrows indicate the shift of Fur-bound probes. The apparent binding affinities 
(Kq) were calculated using GraphPad Prism and were represented as mean ± standard error. The gene designations of N. gonorrhoeae F62 were 
assigned according to their homologues in N. gonorrhoeae FA1090. 
doi:1 0.1 371 /journal.pone.0096832.g002 



According to the above analytical criteria, a large number of 
tested predicted Fur boxes (14 out of 18) in the promoter regions of 
iron-induced gonococcal genes displayed specific binding proper- 
ties to the Fur protein, while 4 predicted Fur boxes failed to bind 
to Fur. SpecificaUy, NG01683 and NG01745 did not show 
increased DSH for any concentration of Fur protein (Figure lA). 
The DSHs of predicted Fur boxes in NGO302 and NGO0377 
promoter regions were only ~1 nm (1.0±0.06 nm and 
1.0 ±0.02 nm, respectively) when 800 nM Fur protein was added 
(Figure lA), suggesting that these two predicted sequences did not 
bind to Fur according to the defined criteria, hi contrast, the other 
14 predicted Fur boxes showed variable binding properties to Fur. 
The predicted Fur boxes in the NGO0436 and NG01738 
promoter region displayed DSH values of 1.4±0.05 nm and 
2.0±0.2 nm, respectively when 800 nM Fur was added 
(Figure lA), corresponding to 0.8 dimer/dsDNA molecule and 
1.6 Fur dimers/dsDNA molecule, respectively (Figure IB). When 
400 nM Fur was added, the predicted Fur boxes in the promoter 
regions of NG01284, NG01845 and NG01948 showed DSH 
values of 1.7±0.03 nm, 1.7±0.3 nm and 1.7±0.07 nm, respec- 
tively, which all corresponded to 0.5 Fur dimers/dsDNA molecule 
(half amount of dsDNA is bound by Fur dimers) (Figure 1 A and 
B). The Fur boxes in the promoter regions of NGO0155, 
NGO0304, NGO0899 and NG01189 had DSH values of 
1.9±0.1 nm, 2.0±0.1 nm, 2.2±0.1 nm and 2.0±0.03 nm, re- 
spectively (Figure lA). These DSH values corresponded to 
approximately one Fur dimer/dsDNA molecule (Figure IB). 



When 200 nM Fur was added, the predicted Fur boxes in 
NGOOlOl, NGO0641 and NG01419 promoter regions had 
DSH values of 1.9±0.4 nm, 2.0±0.1 nm and 1.9±0.4 nm, 
respectively which correlated to one Fur dimer/ dsDNA 
(Figure lA and B). The two highest binding affmities were 
observed for predicted Fur boxes in NGO0073 and NG01957 
promoter regions. When 200 nM Fur was added, DSH values of 
these two Fur boxes were 3.5±0.4 nm and 2.9±0.1 nm, 
respectively, corresponding to 3 Fur dimers and 1.5 Fur dimers/ 
dsDNA molecule, respectively (Figure lA and B). 

To further vahdate Fur binding properties of the predicted Fur 
boxes characterized with IRIS, we performed traditional EMSA 
on two of the newly identified Fur bound probes, NGO0073 and 
NGOOlOl. Fur binding to the predicted Fur boxes in NGO0073 
and NGOOlOl promoter regions was specific as determined by a 
cold competition assay (Figure 2A). Consistent with the binding 
affmity trend estimated using IRIS, the apparent Kjj's of Fur 
binding to Fur boxes of>r gene, NGO0073 and NGOOlOl 
calculated using EMSA results were 31.8±7.9 nM, 39.8±6.3 nM 
and 406.1 ±144.8 nM, respectively (Figure 2B). The combina- 
tion of these results supports IRIS as an accurate method to 
determine the specificity and affinity of Fur binding to an array of 
immobilized dsDNA sequences. 

Of interest among the 18 tested predicted Fur boxes in the 
promoter regions of iron-induced genes, fourteen were shown to 
bind to Fur specifically in vitro (Table 1). No information about 
Fur binding to these promoter regions was known previously. 



PLOS ONE I www.plosone.org 



5 



May 2014 | Volume 9 | Issue 5 | e96832 



A High-Throughput Method for Detecting Protein-DNA Interactions 



NGO0073 







25 ■ 




U. 




+ 




CD 
> 


> 






o 


15 ■ 


< 






RN 




10 - 


E 


o 






g 


5 - 










cr: 


0 ■ 



NA 


7a 1 


en 


■4—' 


E 


o 




gO.5 




iff 

(0 




Of 




WT fiir mutant ftir comp 
NGO0101 



WT ft/r mutant fiircomp 
NGO0155 



WT ftir mutant fiircxjmp 

NGO0436 



M 



i 



* 




WT fur mutant ft/rcomp 

NGO0641 



i 



>K 

* 



1—1 * 



1,2 



WT for mutant furcomp 

NG01738 



WT fur mutant furcomp 

NG01957 



0.9 • 



0.6 ■ 



0.3 ■ 





WT ft/rmutant ftircomp 



WT fur mutant ft/rcomp 



WT fiir mutant ftircomp 



Figure 3. Transcriptional regulation patterns of genes determined by quantitative real-time PCR. The RNA samples were purified from 
cultures of the wild-type (WT) and fur mutant and fur complemented strains under Iron-replete (4-Fe, grey bars) or Iron-deplete (-Fe, white bars) 
conditions 1 h after addition of 100 iiM Iron or 150 \iM desferal. The mRNA levels of fbpA and norB, genes that were repressed and activated by Iron- 
bound Fur respectively, were used as controls for Iron and Fur regulation In N. gonorrhoeae. The mRNA levels observed for the five conditions (WT 
strain under -Fe conditions, fur mutant strain under +Fe and — Fe conditions, and fur complemented strain under 4-Fe and — Fe conditions) were 
compared to the value of WT strain under 4-Fe conditions. The final results were represented as mean ± standard deviation. A * Indicates significantly 
different compared to the mRNA level of WT+Fe. A ** Indicates significantly different compared to the mRNA level of WT-Fe. The gene designations 
of W. gonorrhoeae F62 were assigned according to their homoiogues in N. gonorrhoeae FA1090. 
dol:1 0.1 371 /journal.pone.0096832.g003 



highlighting the power of IRIS in identifying novel promoters 
bound by Fur. 

New targets directly regulated by Fur 

We next examined the regulatory role of Fur in transcription of 
iron-induced genes containing a functional Fur box in the 
promoter regions using quantitative RT-PCR of jV. gonorrhoeae 
wild type, /ar mutant, and fur complemented .strains. This analysis 
revealed 5 new Fur regulated genes [NGO0155, NGO0436, 
NGO0641, NG01738 and NG01957] (Figure 3, Table 1). An 
iron bound Fur-repressed gene, JbpA [6,53], was used as a control 
for Fur and iron regulation. We observed a statistically significant 
increase (5-fold) in fhpA transcript levels in the fur mutant strain 
grown under iron-replete conditions as compared to the wild type 
strain grown under iron-replete conditions. AS expected we 



observed restoration of Fur mediated repression in the Jiir 
complemented strain (Figure 3). Transcription of NGO0155 
was up-regulated in the fur mutant strain compared to the wild 
type strain under iron-replete (H-Fe) conditions, and was partially 
restored in the Jur complemented strain (Figure 3), suggesting that 
iron-bound Fur directly represses transcription of NGO0155. 
Transcriptional levels of NGO0436 and NG01738 were reduced 
in the fur mutant strain compared to those in the wild type strain 
and partially restored in the fur complemented strain grown under 
iron-replete conditions (Figure 3), suggesting that NGO0436 and 
NGOl 738 are directly activated by iron-bound Fur. Transcription 
levels of NGO0641 and NG01957 were up-regulated in the Jur 
mutant strain compared to the wild type strain under iron-deplete 
conditions (-Fe) and were restored in the Jur complemented strain 
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Figure 4. Schematic depicting Fur mediated control mechanisms in N. gonorrhoeae as revealed by IRIS. Solid lines with arrowheads 
indicate direct activation of transcription via Fur binding to promoter regions of gonococcal genes. Solid lines with bars indicate direct repression of 
transcription via binding of Fur to promoter regions of gonococcal genes. Ellipsoids indicate transcriptional regulatory proteins other than Fur. 
Transcriptional repression of NGO0641, encoding a methyltransferase (triangle), results in alteration of DNA methylation (*) and subsequently alters 
transcription of a subset of genes. Transcriptional repression of NGO0155, encoding a putative transcriptional regulator, results in alteration of 
transcription of NGO01 55 target genes. A subset of genes contain a Fur box in their putative promoter regions but were not regulated by Fur under 
the used growth conditions in this study are termed silent Fur binding. 
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(Figure 3). This transcriptional pattern was similar to that 
previously observed in apo-¥\ir mediated repression [15,16]. 

Transcriptional levels of nine genes, NGO0073, NGOOlOl 
(Figure 3), NGO0304, NGO0899, NG01189, NG01284, 
NG01419, NG01845 and NG01948, showed no diflFerence in 
expression among the three strains (Table 1, Figure S2). This 
suggests that these genes are not regulated by Fur under the 
growth conditions in this study. However, due to the fact that Fur 
specifically binds to promoter regions of these genes with high 
affinities, we termed this phenomenon "silent Fur binding" in jV. 
gonorrhoeae. 

Discussion 

IRIS provides a simple microarray screening technique for fast 
identification of in vitro protein-nucleotide interactions, which 
substitutes for onerous gel shift assay. This method is able to 
circumvent complicated and tedious operations such as DNA 
manipulation that may be hindered due to limited availability of 
genetic tools for some organisms and can be a complementary 
method for ChlP-seq [54]. Further application of this microarray 
screening method includes optimization of biochemical conditions 
for protein-nucleotide binding, determination of essential nucleo- 
tides in the conserved binding sequences, as well as determination 
of the dissociation constant of a protein-nucleotide interaction 
through the development of a real-time/ dynamic assay. 

In this study we used Fur, a global transcriptional protein in jV. 
gonorrhoeae to demonstrate the application of IRIS in determining 
the protein-nucleotide interactions. As a successful application of 
IRIS in this study, we identified 14 new Fur binding sequences out 
of 18 tested predicted Fur boxes in the intergenic regions within 
the gonococcal genome and 5 of these 1 4 genes were shown to be 
regulated by Fur under the experimental conditions used 
(Figure 3). NGO0155 encodes a small protein of 4.7 kDa that 
contains a DNA binding domain. NGO0436 and NG01738 
encode a putative hydroxymethyltransferase and a NADH 
dehydrogenase I subvmit M, respectively, which may be important 
for the metabolism of jV. gonorrhoeae. NG01957 encodes a putative 
export protein and NGO0641 encodes a DNA methyltransferase 



(ModA13) of a Type III R-M system. MoclA13 methyltransferase 
recognizes 5'-AGAAA-3' and methylates the third adenine [55]. 
Deletion of NGO0641 or phase OFF of NGO0641 in N. gonorrhoeae 
resulted in transcriptional alteration of 17 genes, which suggests 
that the methyltransferase NGO0641 is involved in epigenetic 
regulation in M. gonorrhoeae [55]. While the relationship between 
Fur and phase variation needs further investigation, it is plausible 
that transcriptional repression of NGO064 1 by Fur may change 
the methylation profile of the gonococcal genome and results in 
epigenetic regulation of the gonococcal transcriptome. 

We also hypothesize a new function for the Fur protein (silent 
Fur binding), based on the observation that ~64% of gonococcal 
genes containing a functional Fur box in promoter regions were 
not regulated by Fur under the used experimental conditions 
(Figure 4). This phenomenon of silent Fur binding has been 
recently demonstrated with a novel gonococcal phage repressor 
protein, Npr [48] . Under //; vivo conditions, some Fur binding sites 
may not be accessible due to the binding of other transcriptional 
regulators and nucleoid-associated proteins including H-NS 
(histone-like nucleoid structuring protein), IHF (integration host 
factor), HU (heat unstable protein) and Fis (factor for inversion 
stimulation) [56] . Thus it is possible that a number of gonococcal 
genes that contain a Fur box in the promoter regions are not 
primarily regulated by Fur, but regulated by other transcriptional 
regulators that may respond to different stimuli. In this scenario, 
cross-talk of Fur with other regulatory proteins would result in co- 
regulation of the same targets in response to multiple environ- 
mental signals encountered during gonococcal infection. Such 
condition-specific expression phenomenon was also discovered in 
the FNR regialon in E. coUv'vSl a genome-scale ChlP-seq study [54]. 

Previous studies on Fur regulated genes have been designed to 
identify Fur regulated genes first in a. fur mutant strain followed by 
in vitro Fur binding assays to confirm the direct regulatory role of 
Fur in gene transcription [8,24,57,58,59]. By design, this type of 
analysis identifies genes that are either direcdy or indirectiy 
regulated by Fur. In contrast, IRIS allows for the rapid 
identification of Fur binding DNA sequences in vitro first, followed 
by subsequent transcriptional analyses. IRIS microarray screening 
of protein-nucleotide interaction will be a powerful high-through- 
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put tool to complement transcriptomt- analysis for the identifica- 
tion of regulatory networks of a broad array of bacterial 
transcriptional proteins. 

Supporting Information 

Figure SI Mass density image of microarray surface 
produced using IRIS. (Left) Following incubation with Fur 
protein (800 nM), an image of the oligonucleotide array is 
produced quantifying the surface mass density across the entire 
surface. By comparing the post-incubation densities to those of the 
pre-incubation image, mass changes attributed to Fur binding can 
be easily determined on a spot by spot basis. (Right) The 
differential spot heights (DSHs) were determined from the mean of 
between 500 and 800 total pixels, depending on spot size, used for 
comparison of the circular spot region (green area) and an outer 
background annulus region (red region). Pixels in the image which 
returned large residtial error during the data fitting process (ex: salt 
residue left in the center of each spot - shown here as missing 
pixels within the green spot region) were eliminated using a 
threshold value. 
(TIF) 

Figure S2 Transcriptional regulation patterns of genes 
determined by quantitative real-time PGR. The RNA 

samples were purified from cultures of the wild-type (WT), Jur 
mutant dtnd Jur complemented strains under iron-replete (H-Fe, grey 
bars) or iron-deplete (-Fe, white bars) conditions 1 h after addition 
of 100 )j.M iron or 150 |J.M desferal. The mRNA levels observed 
for the five conditions (WT strain under — Fe conditions, fur 
mutant strain under H-Fe and — Fe conditions, and Jur comple- 
mented strain under +Fe and — Fe conditions) were compared to 
the value of WT strain under -l-Fe conditions. The final results 



were represented as mean ± standard deviation. A * indicates 
significantly difiFerent compared to the mRNA level of WT-nFe. 
The gene designations of N. gonorrhoeae F62 were assigned 
according to their homologues in jV. gonorrhoeae FA1090. 
(TIF) 

Table SI Experimentally determined Neisseria! Fur 
boxes used as templates for prediction of Fur boxes in 
the genome of N. gonorrhoeae. 

(DOCX) 

Table S2 Predicted Fur boxes in the promoter regions 
at N. gonorrhoeae iron-repressed genes. 

(DOCX) 

Table S3 Double strand DNA probes used in IRIS. 

(DOCX) 

Table S4 Primers used for quantitative RT-PCR. 

(DOCX) 

Table S5 Predicted Fur boxes in the genome of N. 
gonorrhoeae. 

(DOCX) 
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